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Spatial  variations  in  small  scale  vertical  shear  in  the  upper  ocean 
are  described,  relationships  between  small  scale  vertical  shear  and  density 
Stratification  are  investigated,  and  the  potential  for  predicting  mean  vertical 
shear  from  measurements  of  the  density  field  is  assessed.  Fifteen  vertical 
profiles  of  density  and  horizontal  velocity  are  used.  These  were  obtained 
in  several  dynamically  different  ocean  areas  with  a  sensor  package  known 
as  YVETTE.  It  is  concluded  that  a  reasonable  first  approximation  to  a 
climatology  of  small  scale  shear  (at  least,  in  the  North  Atlantic)  might _ 
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ABSTRACT 


Spatial  tanations  in  small  scale  vertical  shear 
in  the  upper  ocean  are  described,  relationships  between 
small  scale  vertical  shear  and  density  stratification  are 
investigated,  and  the  potential  for  predicting  mean  vertical 
shear  from  measurements  of  the  density  field  is  assessed. 
Fifteen  vertical  profiles  of  density  and  horizontal  velocity 
are  used.  These  were  obtained  in  several  dynamically 
different  ocean  areas  with  a  sensor  package  known  as  YVETTE. 
It  is  concluded  that  a  reasonable  first  approximation  to  a 
climatology  of  small  scale  shear  (at  least,  in  the  North 
Atlantic)  might  be  derived  from  a  characterization  of  the 
large  scale  density  stratification.-  Contributions  by 
factors  other  than  density  stratification  to  the  observed 
levels  of  small  scale  vertical  shear  in  the  upper  ocean 


are  also  considered. 


Vie  recommend  that  additional  measure¬ 


ments  be  made  with  an  YVETTE-type  profiler  in  a  wide 
variety  of  ocean  areas.  This  work  was  supported  through 
contract  N0001 4-79-C-0652  from  the  Naval  Oceanographic 
Research  and  Development  Activity. 
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Section  1 
INTRODUCTION 


1.1  BACKGROUND 

Vertical  shear  is  a  parameter  derived  from  the 
vertical  distribution  of  horizontal  velocity,  and  its  impor¬ 
tance  to  Navy  programs  has  been  described  in  a  number  of 
documents.  These  are  cited  and  summarized  in  an  SAI  report 
to  NORDA  by  Lambert  and  Patterson  (1980),  in  which  the 
current  status  of  measurements  of  vertical  shear  in  the 
upper  ocean  was  reviewed.  As  a  result  of  that  report,  it 
was  decided  to  focus  attention  on  shear  over  the  scale  range 
between  1  meter  and  a  few  tens  of  meters.  Two  free-fall 
profiling  instruments  capable  of  resolving  the  velocity 
field  in  this  range,  were  identified  in  that  report.  These 
are  YVETTE  and  SCIMP.  Two  subsequent  SAI  reports  to  NORDA 
have  compared  the  response  of  these  two  instrument  systems 
(Rubenstein  et  al .  ,  1980)  and  presented  profiles  of  data 
obtained  by  them  (Lambert  a_l.,  1980).  Because  of  the 
availability  of  the  YVETTE  data  and  the  wide  variety  of 
ocean  areas  in  which  it  had  been  deployed,  it  was  decided  to 
focus  on  the  YVETTE  data  in  order  to  investigate  the  spatial 
and  oceangrapluc  variability  c upper  ocean  vertical  shear. 

1.2  OBJECTIVES 

This  report  focuses  on  fifteen  vertical  profiles 
of  density  and  horizontal  velocity  obtained  with  YVETTE  in 
several  dynamically  different  ocean  areas.  The  objectives 
of  the  report  are: 
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•  to  describe  the  observed  spatial  variations  in 
small  scale  vertical  shear  in  the  upper  ocean, 
and 


•  to  investigate  the  possible  occurrence  o  f 
systematic  relationships  between  small  scale 
vertical  shear  and  density  stratification. 

1.3  APPROACH 

The  first  step  is  to  become  familiar  with  the 
vertical  distribution  of  the  various  observed  and  derived 
parameters.  Section  2  contains  a  brief  discussion  of  the 
profiles  from  a  single  sample  station,  while  profiles  from 
the  other  stations  are  presented  in  Appendix  A.  It  is  shown 
that  the  magnitudes  of  gradient  quantities  (particularly 
vertical  shear)  in  the  ocean  are  naturally  arranged  into 
distinct  stratification  regimes.  Profiles  from  dynamically 
different  ocean  areas  are  compared  and  contrasted  in  terms 
of  these  stratification  regimes. 

In  Section  3  the  magnitudes  of  gradient  quantities 
within  each  stratification  regime  of  the  sample  station  are 
examined  with  the  use  of  histograms.  Histograms  for 
the  remaining  stations  are  presented  in  Appendix  B.  The 
distributions  revealed  by  the  histograms  for  the  different 
stratification  regimes  are  contrasted  and  discussed.  Since 
the  gradient  quantities  are  computed  by  taking  simple- 
vertical  differences,  the  effects  of  changing  the  vertical 
differencing  interval  are  illustrated  and  explained.  The 
histograms  from  various  ocean  areas  are  compared  and  the 
variation  of  mean  Brunt-Vaisala  frequency  squared  (N2)  and 
mean  vertical  shear  squared  (S^)'  from  one  ocean  area  to 
another  is  illustrated  and  discussed. 
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Relationships  between  N2  and  S2  are  examined 
m  Section  4.  This  includes  an  examination  of  the  effects 
of  varying  the  vertical  differencing  interval  and  the 
vertical  smoothing  interval  ^filter  length)  on  the  correla¬ 
tion  between  N2  and  S2.  The  potential  for  predicting 
S2  from  measurements  of  N2  is  discussed. 

1.4  YVETTE 

At  the  time  these  profiles  were  obtained 
(1973-79).  YVETTE  consisted  of  a  four-meter-long  tube 
equipped  with  Neil  Brown  conductivity,,  temperature,  and 
depth  (CTD)  sensors  as  well  as  an  orthogonal  pair  of  acous¬ 
tic  current  meters  (Figure  1.1)*  A  detailed  descrintion  of 
YVETTE  is  provided  ny  Evans  et  &1.  (1979).  When  deployed,, 
the  probe  falls  freely  through  the  water  column  at  a  rate  of 
20  to  25  cm  s-3.  The  water  temperature,  conductivity,  two 
orthogonal  components  of  horizontal  velocity  relative  to  the 
probe  housing,,  and  the  orientation  of  the  instrument  rela¬ 
tive  to  the  earth's  magnetic  field  are  internally  recorded 
in  digital  form  on  magnetic  cassette  tape.  The  environment 
is  sampled  at  a  rate  of  2,5  times  per  second  to  yield  a 
vertical  sampling  interval  of  8  to  10  cm.  The  CTD  data  have 
a  least  bit  resolution  of  0.001  mmho  cm-3-,  0.0005°C,  and 
0.05  dbar,  respectively.  The  acoustic  current  meters  have  a 
least  bit  resolution  in  speed  of  0.05  cm  s-3  and  the 
compass  can  resolve  orientation  to  within  3°.  Noise  levels 
in  all  sensors  are  dominated  by  least  bit  resolution. 

The  data  set  provided  to  us  by  David  Evans  (Uni¬ 
versity  of  Rhode  Island)  has  been  filtered  to  remove  varia¬ 
tions  on  vertical  scales  less  than  2m.  The  velocity 
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profiles  have  also  been  corrected  for  horizontal  drag 
effects  imposed  on  the  housing  of  YVETTE  by  the  large  scale 
vertical  shear.  A  description  of  this  processing  is  con¬ 
tained  tn  the  SAI  report  to  NORDA  by  Lambert  <at  al_.-  (1980). 
The  resulting  data  set  can  resolve  horizontal  velocity 
structure  with  vertical  scales  between  2m  and  100m.  The 
Instrument  response  to  larger  scale  shear  cannot  be 
estimated. 


1.5  DISTRIBUTION  OF  PROFILES 

The  position,  date,  and  other  relevant  information 
about  the  individual  YVETTE  stations  are  presented  in 
Table  1.1.  The  various  ocean  areas  sampled  by  YVETTF  are 
also  indicated  on  the  map  shown  in  Figure  1.2. 
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Table  1.1 
YVETTE  STATIONS1 


Station 

Number 

Time 

(GMT) 

Date 

Latitude 

GO 

Longitude 

(W) 

Garment 

5 

1942 

5  Nov.  75 

32°19' 

64 “34' 

Near  Bermuda 

8 

0225 

8  Nov.  75 

35*00' 

66 "30' 

Sargasso  Sea 

9 

1219 

It 

" 

" 

it 

10 

1814 

9  Nov.  75 

3S°09' 

69*06' 

Gulf  Stream 

11 

0036 

10  Nov.  75 

38°05' 

69*03' 

li 

12 

1312 

II 

38*15' 

69*07' 

" 

IS 

— 

7  May  77 

22*47’ 

70*43' 

Edge  of  tbemvocline  eddy 

21 

— 

9  May  77 

22*27' 

70°57' 

Center  of  thermocline  eddy 

23 

— 

16  May  77 

36 “24’ 

67*36’ 

Cuter  part  of  GSR2 

24 

— 

•' 

36*20' 

67*44’ 

Midway  along  radius  of  GSR 

25 

— 

17  May  77 

3G°09' 

67*53' 

Near  center  of  GSR 

NOR  1 

-- 

1973 

— 

-- 

Norwegian  Fjord 

NOR  4 

— 

" 

— 

— 

" 

NOR  6 

— 

" 

— 

— 

Norwegian  Coastal  Current 

EP0C336 

— 

July  79 

0*03’ 

109*57' 

Data  start  just  below  EDO1 

1  Adapted  from  Lambert  et  aJL.  (1980) 

2  Gulf  Stream  Ring 

2  Equatorial  Pacific  Ocean  Climate  Study 
4  Equatorial  Undercurrent 
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Section  2 
VERTICAL  PROFILES 

2.1  TYPES  OF  PROFILES 

The  types  of  vertical  profiles  constructed  for 
each  YVETTE  station  are  illustrated  in  Figures  2.1,  2.2,  and 
2.3,  using  YVETTE  station  8  as  an  example.  These  figures 
exhibit  various  aspects  of  the  vertical  distribution  of 
density,  horizontal  velocity,  and  Richardson  number, 
respectively.  Profiles  from  the  other  YVETTE  stations 
are  shown  in  Appendix  A.  The  presentation  here  differs 
from  that  in  Lambert  et^  jU .  (  1980)  in  that  all  profiles 
are  presented  on  a  uniform  set  of  coordinate  scales  to 
facilitate  comparison  between  stations. 

The  vertical  distribution  of  temperature,  salin¬ 
ity,  density  (sigma-t),  Br unt- V a i sal  a  (B.V.)  frequency, 
and  the  vertical  gradient  of  temperature  are  shown  in 
Figure  2.1,  Sigma-t  (jt)  is  related  to  density  (;)  by 

;t  «  (p-1)  *  io3  . 

B.V.  frequency  (N)  is  related  to  the  vertical  gradient  of 
density  by 


where  g  is  gravitational  acceleration  and  z  increases  down¬ 
ward.  The  vertical  gradients  were  computed  by  taking  simple 
differences  along  the  filtered  profiles  over  a  vertical 
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Figure  2.1  Profiles  of  (a)  temperature,  fb)  salinity, 

(c)  density  (ct)>  (d)  Brunt -Vais ala  frequency, 
and  (e)  vertical  gradient  of  temperature  for 
YVETTE  station  8  occupied  on  8  November  1975 
at  35"00'N,  66°30'W  (Sargasso  Sea).  The  pro¬ 
files  can  be  divided  into  stratification  regimes 
as  indicated  by  the  dotted  lines  at  135  and  250  dbars. 
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Figure  2.2  Profiles  of  the  (a)  east  and  (b)  north  component  of 
velocity,  the  (c)  ..agnitude  and  '(d)  direction  of  the 
horizontal  velocity  vector,  and  the  vertical  gradient 
of  the  (e)  east  and  (f)  north  component  of  velocity 
for  YVETTE  station  8  m  the  Sargasso  Sea.  The  dotted 
lines  are  at  135  and  250  dbars. 


(a) 


(c) 


Figure  2.3  Profiles  of  (a)  Brint-Vaisal'a  frequency  squared, 

N2 ,  (b)  vertical  shear  squared,  S2,  and  (c)  Logjo 
(Ri),  where  Ri  =  N2/S2,  for  YVETTE  station  8  in 
the  Sargasso  Sea  The  vertical  line  in  (c)  is  at 
Ri  =  0.25.  The  dotted  lines  are  at  135  and  250  dbars. 


interval  of  approximately  2  dbars.  For  the  purposes  of  this 
paper,  1  dbar  is  considered  to  be  equal  to  1  meter..  In 
regions  of  near  neutral  stability,  such  as  within  the 
surface  mixed  layer,  noise  in  the  temperature  and/or  salin¬ 
ity  profile  can  produce  apparent  inversions  in  the  density 
profile.  These  apparent  inversions  result  in  imaginary  B.V. 
frequencies.  As  a  matter  of  convenience  in  constructing  the 
B.V.-  profiles,  the  imaginary  B.V.  frequencies  have  been 
plotted  as  negative  real  quantities. 

The  water  column  sampled  at  station  8  can  be 
divided  into  three  stratification  regimes  on  the  basis  of 
the  vertical  gradient  of  either  temperature  or  density. 
This  division  is  indicated  in  Figures  2.1,  2.2,  and  2.3  by 
dashed  lines  at  135  dbars  and  250  dbars.  The  placement  of 
the  boundaries  between  regimes  is  somewhat  subjective. 
Above  135  dbars  is  the  surface  mixed  layer  (ML),  character¬ 
ized  by  nearly  constant  temperature,  salinity,  and  density. 
Between  135  dbars  and  250  dbars  is  the  near-surface  thermo- 
cline  (ST)  characterized  by  a  relatively  large  mean  vertical 
gradient  in  temperature  and  density.  Below  the  thermocline 
is  a  regime  (henceforth  abbreviated  as  BST)  of  significantly 
smaller  mean  vertical  gradients  in  temperature  and  density. 
We  note  that  in  the  case  of  YVETTE  8,  the  BST  regime  could 
be  further  divided.  Between  approximately  300  and  400m 
there  is  a  nearly  homogeneous  layer  (18°C  water)  which  is 
characterized  by  uniformly  low  B.V.  frequency.  Below  this 
layer  is  the  top  of  the  main  thermocline  where  the  tempera¬ 
ture  begins  to  decrease  and  the  B.V.  frequency  begins  to 
increase.  The  profiles  from  the  other  YVETTE  stations  can 
be  similarly  divided  into  stratification  regimes  on  the 
basis  of  the  relative  magnitudes  of  vertical  gradients, 
however,  the  depth  and  thickness  of  the  regimes  vary  from 
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station  to  station  (see  Appendix  A).  Since  the  profile  data 
begin  at  various  distances  below  the  sea  surface,  in  many  of 
the  profiles  no  surface  mixed  layer  is  observed. 

Various  aspects  of  the  vertical  distribution  of 
horizontal  velocity  are  shown  in  Figure  2.2.  Profiles  of 
the  east  (u)  and  north  (v)  components  of  horizontal  velocity 
are  shown  in  panels  (a)  and  (b),  respectively.  Profiles  of 
the  magnitude  (speed)  and  direction  of  the  horizontal 
velocity  vector  are  shown  in  panels  (c)  and  (d).  In  the 
direction  plot,  north  corresponds  to  0°,  east  corresponds 
to  90°,  etc.  The  vertical  gradient  of  the  east  and  north 
components  of  velocity  are  shown  in  panels  (e)  and  (f).  It 
is  evident  that  these  velocity  gradients  exhibit  a  broader 
range  of  values  within  the  ST  regime  than  elsewhere.  This 
means  that  there  are  sharper  or  more  abrupt  variations 
(along  the  vertical  axis)  of  the  horizontal  velocity  vector 
within  the  thermocline  than  elsewhere.  These  variations  can 
be  due  to  changes  in  either  the  magnitude  or  the  direction 
of  the  '  tor,  or  both. 

The  variations  in  du/dz  and  dv/dz  are  combined  and 
made  non-negative  by  computing  shear  squared  (S2),  where 


The  gross  similarity  between  the  profile  of  B.V.  frequency 
squared  (N2)  and  that  of  S2  is  seen  by  comparing  panels 


» 
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(a)  and  (b)  of  Figure  2.3.  Both  profiles  exhibit  rela¬ 
tively  small  values  within  the  ML  and  BST  regimes,  but 
relatively  large  peaks  within  the  ST  regime.  This  suggests 
that  some  relationship  may  exist  between  N2  and  S2 . 
However,  the  profile  of  the  logarithm  of  Richardson  number 
(Ri),,  where  Ri  =  N2/S2,  exhibits  a  broad  range  of  values 
in  panel  (c)  of  Figure  2.3.  Hence  a  linear  relationship 
between  N2  and  S2  that  is  suggested  by  comparing  whole 
profiles  (that  is,  500  m  profiles}  does  not  hold  when 
comparing  individual  points.  The  vertical  line  in  panel  (c) 
at  Logio  (Ri)  =  -0.602  corresponds  to  Ri  =  0.25.  Simple 
linear  theory  predicts  that  Ri  must  be  less  than  this 
critical  value  for  the  existence  of  shear-induced  instabil¬ 
ities.  Figure  2.3  suggests  that  most  of  the  very  low  values 
of  Ri  occur  in  portions  of  the  water  column  where  N2  is 
low  rather  than  in  portions  where  S2  is  high.  This 
same  tendency  is  evident  in  all  YVETTE  profiles  (see 
Appendix  A). 

In  regions  of  very  low  S’ 2  or  S2 ,  the  noise  level 
inherent  with  the  sensor  package  causes  the  relative  error 
in  the  calculated  values  of  N2  and  S2  to  become  large. 
Therefore,,  in  such  regions  the  relative  error  in  Ri  is 
also  quite  large.  From  the  precision  of  the  various  sensors 
we  have  estimated  (Molinelu  e£  al_.  ,  1981)  that  the  signal- 
to-noise  ratio  in  the  N2  profile  is  greater  than  unity 
when  N2  >  1.8  x  10-8  s-2,  and  that  ir,  the  S2  profile 
is  greater  than  unity  when  S2  >  3.3  x  10-8  s-2.  Outside 
the  surface  mixed  layer,  the  values  of  N2  and  S2  are  only 
rarely  below  the  estimated  noise  level.  A  more  detailed 
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discussion  of  the  effects  of  noise  leve.  ;n  YVETTE  data  on 
calculations  of  Ri  is  presented  in  Newman  et  al ■  (1981). 
While  YVETTE  is  descending  through  the  water  column,  various 
dynamic  responses  of  the  system  may  increase  the  .'vise  level 
above  that  estimated  here.  This  issue  remains  to  be 
resolved.  In  any  case.-  the  actual  noise  level  is  nut 
expected  to  have  any  appreciable  impact  on  the  gross 
characterizations  presented  herein. 

2.2  COMPARISONS  BETWEEN  STATIONS 


Representative  profiles  from  each  ocean  area 
sampled  by  YVETTE  are  presented  in  Figure  2.4  to  contrast 
the  different  oceanic  conditions  encountered.  It  can  be 
seen  that  a  fairly  broad  range  of  oceanic  conditions,  from 
the  cold,  fresh  waters  off  the  coast  of  Norway  to  the  warm, 
salty  waters  of  the  Sargasso  Sea,  are  represented  by  the 
data  set.  Waters  in  the  equatorial  Pacific  are  colder  and 
fresher  than  most  of  the  waters  of  the  North  Atlantic. 
Although  the  temperature  and  salinity  profile  from  the 
Gulf  Stream  (station  10)  exhibits  substantial  interleaving, 
the  sigma-t  profile  shows  that  these  fluctuations  are 
essentially  density-compensating.  The  amplitude  of  the 
temperature  and  salinity  interleaving  increases  from  south¬ 
east  to  northwest  across  the  Gulf  Stream.  (Compare  stations 
11,  10,  and  12,  in  sequence,  in  Appendix  A.)  We  note 
from  the  sigma-t  profiles  in  Figure  2.4  that  below  50  meters 
density  increases  along  constant  pressure  surfaces  radially 
outward  from  a  point  in  the  central  Sargasso  Sea  (station  8) 
toward  Bermuda  (station  5),  the  Bahamas  (station  18),  and 
the  Gulf  Stream  (station  10).  This  is  due  to  the  large 
scale  anticyclomc  circulation  in  the  Sargasso  Sea. 
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A  profile  of  shear  squared  for  YVETTE  station  8 
is  shown  in  Figure  2.3.  Profiles  of  shear  squared  for  the 
other  representative  stations  are  shown  in  Figure  2.5.  It 
is  evident  that  the  vertical  distribution  of  shear  squared  is 
just  as  diverse  as  any  of  the  other  parameters. 

The  contrasts  from  one  stratification  regime 
to  another  and  from  one  ocean  area  to  another  will  be 
further  illustrated  and  discussed  in  Sections  3  and  4. 
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Section  3 
HISTOGRAMS 


3.1  INTRODUCTION 

The  YVETTE  profiles  presented  in  the  preceeding 
section  showed  the  distribution  of  various  parameters  with 
depth..  It  was  noted  that.-  for  each  station,  the  profiles  of 
the  gradient  quantities  can  be  divided  along  the  depth  axis 
into  two  or  more  distinctly  different  regimes  which  cor¬ 
respond  to  regimes  of  different  mean  density  stratification. 
In  order  to  further  examine  the  differences  between  these 
stratification  regimes  and  the  differences  between  ocean 
areas,  histograms  of  the  gradient  quantities  were  con¬ 
structed  for  each  regime  and  each  different  ocean  area. 

The  data  points  in  the  processed  YVETTE  profiles 
are  equally  spaced  in  time.  Due  to  variations  in  the 
descent  rate  of  the  instrument,  the  points  are  not  always 
equally  spaced  in  depth,  although  the  nominal  spacing  is  lm. 
The  gradient  value  at  each  observed  deptc.  was  computed  by 
differencing  the  values  above  and  below  the  specified  depth 
and  dividing  by  the  differencing  interval,  iz.  Since  the 
gradients  are  computed  in  this  way,  the  effects  of  changing 
the  vertical  differencing  interval  are  illustrated  and 
expl  am'  •* 


3.2  TYPES  OF  HISTOGRAMS 

As  in  Section  2,  YVETTE  8  is  used  to  illustrate 
the  types  of  histograms  constructed  (Figures  3.1,  3.2.  3.4, 
and  3.7).  Histograms  of  B.V.  frequency  squared  ( N  2 )  , 
shear  squared  (S2),  and  Richardson  number  (Ri)  for  the 
other  YVETTE  stations  are  presented  in  Appendix  B. 
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Note  that  the  vertical  axis  of  these  histograms 
is  in  units  of  observed  probability  or  fractional  frequency 
of  occurrence.  When  the  observed  probability  exceeds  0.35, 
the  histogram  bar  is  truncated  and  the  observed  probability 
is  enclosed  in  parentheses  above  the  appropriate  bar.  For 
example,  in  Figure  3.1  (a),  0.7  of  the  observed  values  occur 
in  the  interval  between  0.0  and  0.05  x  10-4  rad2  s-2. 
While  the  observed  range  of  values  for  each  parameter  is 
quite  large,  in  many  of  the  distributions  the  majority  of 
the  values  are  contained  within  a  relatively  narrow  range. 
We  note  that  the  shape  of  a  histogram  can  be  quite  sensitive 
to  bin  width.  Hence,  the  range  and  bin  widths  were  selected 
to  accommodate  most  of  the  values  observed  in  the  entire 
data  set  (i.e.,  15  stations)  and  at  the  same  time  provide  a 
reasonably  smooth  distribution  m  the  vicinity  of  the  mode 
(i.e.  the  peak  or  most  probable  value).  Once  selected,  the 
range  and  bin  width  for  each  parameter  are  held  constant  so 
that  the  histograms  can  be  easily  compared.  Values  that 
exceed  the  bounds  of  the  horizontal  axis  are  accumulated  m 
the  last  bin  on  the  extreme  left  or  right,,  whichever  is 
appropriate.  For  example,  in  Figure  3.1  (b).  0.42  of 
all  values  are  greater  than  1.3  x  10~4  rad2  s-2. 

The  mean,  median,  and  variance  of  each  distribu¬ 
tion  are  indicated  m  the  legend  on  each  histogram.  The 
depth  range  and  number  of  observations  within  each  stratifi¬ 
cation  regime  are  indicated  in  the  legends  of  the  histograms 
on  the  left  side  of  each  figure.  Since  a  gradient  value  is 
computed  at  each  observed  depth,  the  number  of  observations 
is  approximately  one  per  meter.  In  the  histograms  of  Ri 
(e.g.  Figure  3.7)  the  legend  includes  a  line  indicating  the 
observed  probability  of  Ri  less  than  the  critical  Richardson 
number,  0.25. 
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Recall  that  at  station  8  the  water  column 
has  been  divided  into  three  regimes:  the  mixed  layer  (22  to 
135  dbars),  the  near-surface  thermocline  (135  to  250  dbars), 
and  the  region  below  the  thermocline  (250  to  493  dbars)’. 
The  observed  distributions  of  density  gradient  and  tempera¬ 
ture  gradient  values  in  these  three  regimes  are  shown  in 
Figure  3.1.  Panels  (a),  (b),  and  (c)  show  histograms  of 
N2  for  the  ML,  ST,  and  3ST  regimes,  respectively.  (Com¬ 
pare  with  panel  (a)  in  Figure  2.3.)  Panels  (d),  (e),  and 
(f)  show  the  corresponding  histograms  of  dT/dz.  (Compare 
with  panel  (e)  in  Figure  2.1). 

Since  temperature  typically  decreases  with  depth 
and  density  increases  with  depth,  the  majority  of  the 
values  of  dT/dz  are  opposite  in  sign  to  those  of  B.V. 
frequency  squared.  To  the  extent  that  density  variations 
are  dominated  by  variations  in  temperature  rather  than 
salinity,  the  histograms  of  dT/dz  should  be  the  mirror  image 
of  those  for  B.V.  frequency  squared.  The  contrast  between 
the  three  stratification  regimes  is  clearly  illustrated  by 
the  histograms.  The  majority  of  the  values  within  the  ML 
are  near  zero.  Within  the  BST  regime,  the  water  column  is 
more  stably  stratified  and  the  variance  of  the  distribution 
is  somewhat  larger.  Within  the  ST  the  stratification  is 
much  more  stable  and  the  variance  of  the  distribution  is 
quite  large.  This  same  pattern  is  evident  in  the  histograms 
of  N2  for  the  other  YVETTE  stations  shown  m  Appendix  B. 
Hayes  et  al^.  (  1975)  have  noted  that  a  Poisson  distribution 
function  provides  a  reasonably  good  visual  fit  to  a  similar 
set  of  histograms  of  dT/dz. 
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Histograms  of  du/dz  art-  presented  in  Figure  3.2. 
Once  again  the  contrast  between  the  three  stratification 
regimes  is  clearly  illustrated  by  the  histograms.  (Compare 
with  panel  (e)  in  Figure  2.2.)'  The  \ariance  is  signif¬ 
icantly  greater  within  the  ST  than  within  the  ML  or  BST 
regimes.  Panels  (a),  lb),  and  (c)  show  distributions 
of  du/dz  computed  by  differencing  over  a  vertical  interval 
of  2m.  Panels  ( d  5 ,  (e),  and  (f)  were  computed  by  differ¬ 
encing  over  8m.  Increasing  the  vertical  differencing 
interval,  iz,  effectively  subsamples  the  velocity  profile 
provided  by  YVETTE.  This  is  analogous  to  measuring  the 
velocity  profile  with  a  lower  resolution  discrete  tech¬ 
nique  such  as  a  vertical  array  of  current  meters.  Of 
course,  such  discrete  data  are  aliaced  to  some  extent 
by  small  scale  structure.  A  lower  resolution  profiling 
instrument  would  not  be  sensitive  to  the  small  scale  struc¬ 
ture.  Such  profiler  data  could  be  simulated  by  low-pass 
filtering  the  YVETTE  profiles  to  remove  the  high  wavenumber 
structure. 

Although  the  du/dz  histogram  shapes  resemble  a 
Gaussian  probability  distribution  function,  similar  histo¬ 
grams  (Figure  3.3)  statistically  tested  by  Simpson  (1975) 
proved  to  be  nor.-Gaussi an  to  a  high  level  of  confidence.  It 
is  seei  in  Figure  3.2  that  the  shape  of  the  distributions  of 
du/dz  is  sensitive  to  changes  in  the  vertical  differ¬ 
encing  interval.  For  example,  in  all  regimes  the  variance 
of  the  distributions  decreases  as  i z  increases.  This  is 
because  the  contributions  to  the  "tails''  of  the  du/dz 
histograms  come  chiefly  from  the  small  scale  structure 
in  the  velocity  profile.  When  Lz  increases,  the  ability  to 
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Figure  3.3  Histograms  of  the  velocity  shear  components  for  the 
depth  range  66S-S06  m  at  o2  23.2  N,  16  1S.0 

(near  Rockall  Trough).  The  continuous  curve  represents 
the  normal  distribution.  Each  histogram  contains  770 
data  points  and  the  abscissa_is  in  units  of  one  standard 
deviation  (from  Simpson.  1975) 
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resolve  this  structure  is  lost,  and  hence  variance  decreas¬ 
es..  Lower  vertical  resolution  within  a  single  profile  also 
implies  a  smaller  number  of  statistically  independent 
observations,  which  reduces  the  confidence  in  any  statisti¬ 
cal  tests  imposed.  Therefore,  the  details  of  the  distribu¬ 
tion  shape,  such  as  whether  or  not  it  fits  a  Gaussian 
probability  curve,  may  depend  as  much  on  the  vertical 
resolution  of  the  velocity  profile  as  on  the  ocean  environ¬ 
ment  itself.  The  results  of  statistical  tests  of  the  YVETTE 
data  are  presented  by  Hebenstreit  and  Grabowski  (1981). 

Histograms  of  vertical  shear  squared  (S2)  are 
presented  in  Figure  3.4.  As  with  the  du/dz  histograms, 
distributions  of  shear  squared  computed  using  a  vertical 
differencing  interval  of  2m  are  presented  m  panels  (a), 
(b),  and  (c)  and  those  computed  using  an  interval  of  8m  are 
presented  in  panels  (d),  (e),  and  (f).  The  difference 
between  stratification  regimes  is  as  expected  from  panel 
(b)  Figure  2.3.  That  is,  small  values  of  S2  are  more 
common  in  the  ML  and  BST  regimes  than  within  the  ST  regime. 
This  is  a  general  result  as  illustrated  by  the  S2  histo¬ 
grams  (computed  using  a  vertical  differencing  interval  of 
2m  only)  from  the  other  YVETTE  stations  presented  in 
Appendix  B.  The  S2  histograms  resemble  Chi-square  prob¬ 
ability  distnoution  functions.  If  the  shear  components  do 
follow  Gaussian  distributions,  then  this  observation  is  not 
surprising.  (The  distribution  of  the  sum  of  the  squares  of 
two  normally  distributed  random  variables  is  the  Chi-square 
distribution.)  Furthermore,  this  suggests  that  Chi-square 
distributions  might  serve  as  suitable  models  of  S2  distri¬ 
butions.  This  line  of  inquiry  is  addressed  by  Hebenstreit 
and  Grabowski  (1981). 
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Figure  3.4  Histograms  of  vertical  shear  squared  (S2)  computed  by 
differencing  over  a  vertical  interval  of  2  m  for  (a) 
the  ML,  (b)  the  ST,  and  (c)  the  3ST  regime.  Histograms 
of  S2  computed  by  differencing  over  an  interval  of  8  m 
for  (d)  the  ML,  (e)  the  ST,  and  (f)  the  BST  regime. 


Since  the  shape  of  the  du/dz  histogram  is  sen¬ 
sitive  to  the  vertical  differencing  interval,,  the  same  is 
necessarily  true  for  histograms  of  S2.  It  is  seen  in 
Figure  3.4  that  in  all  regimes,  both  the  mean  and  the 
variance  of  the  distributions  decrease  as  Az  increases.  The 
variation  of  mean  S2  with  stratification  regime  and  with  Az 
is  further  illustrated  in  Figure  3.5.  The  effect  of 
increasing  Az  on  the  mean  value  of  S2  distributions  is  a 
general  result  which  is  illustrated  m  Figure  3.6.  This 
effect  has  been  modeled  by  Grabowski  (1980).  These  rela¬ 
tions  suggest  that  it  may  be  possible  to  use  some  of  the 
lower  resolution  velocity  profilers,  such  as  the  EMVP 
(Sanford  et,  al .  ,  1978),  the  XCP  (Sippican  Corp.  ,  1979),  or 
even  vertical  arrays  of  current  meters  (Pollard,  1980),  to 
provide  at  least  rough  estimates  of  the  mean  S2  that  would 
have  been  observed  by  a  high  resolution  profiler  like 
YVETTE.  This  might  allow  data  sets  obtained  with  different 
resolution  instruments  to  be  compared.  A  discussion  of  many 
of  the  existing  velocity  profilers  is  presented  in  Lambert 
and  Patterson  (1980), 

As  noted  earlier  the  ratio  of  N’2  to  S2  is  the 
Richardson  number  (Ri).  Histograms  of  Ri  are  presented 
in  Figure  3.7.  The  format  is  the  same  as  for  Figure  3,2  and 
3.4  except  that  the  observed  probability  of  subcritical  Ri 
(i.e.,  Ri  <  0.25)  is  indicated  m  the  legend.  It  was  noted 
in  Section  2.1  that  when  considering  the  profiles  as  a 
whole,,  N2  and  S2  appear  to  be  related  to  each  other.  How¬ 
ever,  the  profile  of  Log^o  (Ri)  presented  in  Figure  2.3, 
rather  than  being  nearly  constant,  exhibited  a  broad  range 
of  values.  From  examination  of  the  histograms  (Figures  3.1 
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Figure  3.6  Mean  vertical  shear  squared  (S-)  versus 

differencing  interval  (  z)  for  all  YVETTE 
stations  except  2-1  which  is  a  very  short 
profile.  Here  the  mean  in  each  case  is 
computed  using  data  from  the  entire  profile 
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*  Figure  3.7  Histograms  of  Richa-  son  number  (Ri)  computed  by 

differencing  over  a  vertical  interval  of  2  m  for 
(a)  the  ML,  (b)  the  ST,  and  (c)  the  BST  regime. 
Histograms  of  Ri  computed  over  an  interval  of  8  m 
for  (d)  the  ML,  (e)  the  ST,  and  (f)  the  BST  regime. 


and  3.4},  we  once  again  observe  than  N2  and  S2  exhibit 
qualitatively  similar  behavior,  namely,  a  predominance 
of  low  values  within  the  BST  regime  and  a  broad  range  of 
values  within  tne  ST  regime.  In  the  case  of  Ri ,  however, 
the  histrograms  for  the  ST  an'1  R5T  regimes  are  remarkably 
similar.  This  is  as  o 'e  might  expect  from  the  profiles 
shown  in  Figure  2.3.  Henc>,  while  the  distributions  of  N'2 
and  S2  within  the  ST  regime  are  markedly  different  than 
those  below,,  the  distribution  of  their  ratio  remains  rela¬ 
tively  constant.  This  situation  does  not  hold,  however,, 
within  the  ML  regime.  Under  conditions  of  extremely  weak 
density  stratification,  the  presence  of  any  significant 
shear  will  result  in  \ery  low  values  of  Ri .  4s  seen  in 
Figure  3.4,  the  values  of  S2  within  the  ML  are  significant 
and  are  comparable  to  those  within  the  BST  regime,  whereas 
the  values  of  N’2  (Figure  3.1)  in  these  regimes  are 
significantly  different.  Hence,  while  it  seems  quite 
appropriate  on  the  basis  of  N2  or  S2  to  divide  the  water 
column  into  three  regimes,  it  seems  on  the  basis  of  Ri  that 
the  column  should  be  divided  into  only  two  regimes,  namely, 
the  surface  mixed  layer  and  the  region  below  the  mixed 
layer.  This  same  behavior  is  observed  in  the  histograms  for 
the  other  YVETTE  stations  shown  in  Appendix  B. 

Since  histograms  of  S2  are  sensitive  to  changes 
in  the  vertical  differencing  interval  (histograms  of  N2 
are  relatively  insensitive),  it  is  not  surprising  that 
histograms  of  Ri  arc  also  sensitive.  In  Figure  3.7  we 
observe  that  as  Lz  increases  from  2m  to  8m,  the  Richardson 
number  distribution  shifts  to  higher  values.  This  is 
consistent  with  the  shift  of  the  S2  distribution  toward 
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lower  values  (Figure  3.4).-  We  also  note  that  the  change  in 
histogram  shape  induced  by  differencing  over  different 
vertical  intervals  is  substantially  greater  than  the  slight 
contrast  observed  between  the  ST  and  the  BST  regimes.  This 
strongly  suggests  that  the  vertical  resolution  of  the 
profile  of  velocity  is  a  more  important  factor  in  determin¬ 
ing  the  shape  of  the  Ri  distribution  than  the  variations  in 
the  ocean  environment  itself.  This  conclusion  should  be 
borne  in  mind  when  attempting  to  generalize  from  a  single 
data  set  or  when  attempting  to  compare  results  obtained  with 
different  sensors.  Earlier,  we  alluded  to  the  potential  for 
using  values  obtained  with  one  resolution  to  estimate 
the  mean  value  that  would  have  been  observed  with  an 
instrument  with  a  different  resolution.  Such  an  exercise 
might  allow  one  to  compare  Ri  distributions  obtained  with 
different  instruments.  The  effects  of  vertical  resolution 
on  calculated  values  of  Ri  is  illustrated  and  discussed  in 
more  detail  in  Newman  et  al.  (1981). 

3.3  COMPARISONS  BETWEEN  STATIONS 


Histograms  for  YVETTE  8  were  presented  m  the 
preceedmg  section.  Histograms  of  B.V.  frequency  squared, 
vertical  shear  squared,  and  Richardson  number  for  the 
remaining  YVETTE  stations  are  presented  in  Appendix  B.  The 
effect  of  variations  m  the  vertical  differencing  interval, 
iz ,  on  histogram  shape  has  been  illustrated.  Since  this 
effect  is  consistent,  the  histograms  presented  in  Appendix  B 
are  all  constructed  using  a  vertical  differencing  interval 
of  2m.  Clearly,  to  compare  histograms  from  different 
stations,  it  is  necessary  to  hold  ‘.z  constant.  The  differ¬ 
encing  interval  of  2m  was  selected  to  take  advantage  of  fie 
finest  resolution  afforded  by  the  filtered  YVETTE  profiles. 


While  details  of  the  histogram  shapes  (Appendix  B) 
change  from  station  to  station,  the  tendencies  noted  for 
station  8  seem  to  prevail  throughout  the  data  set.  In  some 
'’ases  (e.g.  YVETTE  24,  Figure  B.9)'  the  histogram  shapes  in 
each  stratification  regime  are  poorly  resolved  because  the 
profile  is  short  and  there  are  consequently  only  a  small 
number  of  observations.  At  YVETTE  stations  N 1  and  N4 
(Figures  B.ll  and  B.12)  the  observed  mean  stratification  and 
vertical  shear  squared  within  the  upper  regime  are  so  great 
>5X10-4  s-2)'  that  the  values  are  off  scale.  Surprisingly 
the  Richardson  number  histograms  conform  well  to  those  of 
the  other  stations.  In  one  case,  the  ST  regime  at  YVETTE 
station  E6  (Figure  B.14)',  the  Richardson  number  histogram  is 
somewhat  anomalous.  Here  the  broad  range  of  values  of  B.V. 
fiequency  squared  are  accompanied  by  predominantly  small 
values  of  shear  squared.  This  results  in  unusually  high 
values  (mode  >  2)  of  Richardson  number.  It  is  noted, 
however,  that  this  histogram  is  one  of  those  whose  shape  is 
poorly  resolved  because  of  the  small  number  of  observations. 

Rather  than  comparing  histograms,  the  variations 
from  one  station  to  another  can  be  more  concisely  shown  by  a 
plot  of  a  single  pa-ameter  of  the  histogram,  such  as  the 
mean,  versus  station  '.umber.  Therefore,  the  variations  from 
station  to  station  if  the  mean  value  of  N2  and  S2  (X2 
and  S2)  within  each  stratification  regime  are  shown  in 
Figure  3.8.  In  th  s  figur.  the  stations  are  arranged  in 
descending  order  with  respect  to  the  mean  value  of  S2 
(computed  at  each  s' at  ion  with  values  from  the  whole  pro¬ 
file),.  Since  some  of  the  profiles  are  relatively  short 
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(YVETTE  24,  Nl,  N4,  and  N6),  the  mean  values  for  these 
profiles  are  not  strictly  comparable  to  the  others.  Bearing 
this  in  mind,  they  are  included  anyway  for  the  sake  of 
completeness. 

The  high  mean  values  at  YVETTE  stations  Nl  and  N4 
are  due  to  the  extremely  high  values  in  the  upper  regime 
(i.e.,  the  upper  34m)'  and  are  not  truly  representative  of 
the  bulk  of  the  water  column,  which  is  more  appropriately 
represented  by  the  lower  regime.  (See  profiles  in  Figures 
A. 33  and  A. 36,  respectively.)  Except  for  YVETTE  statons  Nl 
and  N4,  there  is  a  tendency  in  the  North  Atlantic  profiles 
for  S2  levels  to  decrease  with  distance  from  the  Gulf 
Stream,  that  is,  to  decrease  from  the  Gulf  Stream  to  the 
Gulf  Stream  ring,  the  Sargasso  Sea,  Bermuda,  the  Bahamas.- 
and  finally  the  Norwegian  Coastal  Current.  Possibly 
the  same  trend  would  be  observed  in  the  North  Pacific  with 
distance  from  the  Kuroshio.  kith  existing  data,  we  can  only 
observe  that  in  the  equatorial  Pacific,  the  S2  levels  are 
relatively  low. 

Examination  of  Figure  3.8  reveals  that,-  to  a  first 
approximation,  N2  levels  computed  over  whole  profiles 
decrease  with  the  corresponding  S^  levels.  The  anom¬ 
alously  high  N2  levels  at  YVETTE  station  24  are  due  to 
the  very  high  values  within  the  thermocline  which  are 
probably  an  artifact  of  poor  resolution  due  to  the  abrupt 
termination  of  the  profile  (Figure  A. 28).  The  N2  and  S2 
levels  within  different  stratification  regimes  rank  them¬ 
selves  as  expected.  Note  that  the  S2  levels  within 
the  surface  mixed  layer  sometimes  exceed  those  below  the 
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therraocline  (e.g.  YVETTE  stations  8  and  9).  Note  also  that 
within  each  stratification  regime  there  is  a  tendency  for 
the  N2  and  S2  levels  to  vary  up  and  down  in  parallel 
fashion.  This  tendency  will  be  examined  further  in 
Section  4.- 


From  the  temperature  profiles  shown  in  Figures 
A. 22,  A. 25,  and  A. 28,  it  is  evident  that  the  Gulf  Stream 
ring  sampled  by  YVETTE  stations  23,  24,  and  25  is  a  cold 
core,  cyclonic  ring.  While  it  is  not  possible  with  only- 
three  profiles  to  know  the  precise  location  of  the  center 
of  the  ring,  we  note  in  Figure  3.8  the  trend  of  increasing 
N2  and  S2  toward  the  center  of  tms  mesoscale  feature. 
(Recall  that  the  N2  levels  at  YVETTE  24  are  suspect.-) 
The  same  trend  is  evident  in  the  thermocline  eddy  (a  meso¬ 
scale  feature  confined  in  vertical  extent  to  the  thermocline 
portion  of  the  water  column)  observed  near  the  Bahamas  by 
YVETTE  stations  18  and  21.  The  mesoscale  feature  is  most 
clearly  revealed  by  the  lens  of  relatively  high  salinity- 
water  observed  at  station  21  between  100  and  250  dbars.- 
(cf.  Figures  A. 16  and  A. 19.) 

The  tendency  for  N2  to  increase  toward  the 
center  of  a  cyclonic  eddy  may  be  due  to  the  convergence  (in 
the  vertical)  of  surfaces  of  constant  density  toward  the 
center  of  such  features.  This  phenomenon  is  observed  in 
the  vertical  sections  across  Gulf  Stream  rings  shown  in 
Figure  3.9.  (We  note  that  in  the  region  of  the  Gulf  Stream, 
surfaces  of  constant  density  are  almost  parallel  to  surfaces 
of  constant  temperature.)  The  accompanying  increase  in  S2 
toward  the  center  of  such  features  is  not  yet  explained, 
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Figure  3  0  Temperature  sections  between  (a)  Florida  and  3S°N'.  6S°i.' 

(KNORR  65  9-14  Aoril  1077)  ana  (b)  Florida  and  \'»,i  Fn^land 
(ENDEAVOR  11.  26  Julj-18  August  1977)  Several  Gulf  St  rear, 
cjclonic  rings  (as  confirmed  by  snort  X5T  surievs  or  oy  drifter 
trajectories)  were  crossed  including  those  called  Bob,  Charlie, 
and  Dave.  Also  an  anticiclomc  ring  was  observed  northwest  of 
the  Gulf  Ctieam  (from  Richardson,  1980) 
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but.  it  is  consistent  with  the  observed  tendency  for  "s^  to 
vary  with  N^.  From  the  above  discussion  we  would  expect 
both  N**  and  to  decrease  toward  the  center  of  an 
anticyclonic  ring.  Since  we  have  no  data  from  an  anti- 
cyclonic  ring,  this  remains  as  an  hypothesis  to  be  tested 
during  future  field  work. 
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Section  4 


RELATIONSHIPS  BETKEEES  VERTICAL  SHEAR 
AND  DENSITY  STRATIFICATION 

4.1  INTRODUCTION 

In  Figure  2.3,  we  noted  that  a  visual  comparison 
of  the  profiles  of  N2  and  S2  suggests  that  some  rela¬ 
tionship  exists  between  these  two  parameters.  In  this 
figure,  individual  values  were  obtained  by  differencing  over 
a  vertical  interval  of  2m.  However,  individual  values  of 
Ri,  the  ratio  of  \'2  to  S2,  fluctuate  over  several  orders 
of  magnitude,  and  no  relationship  is  evident  when  individual 
values  of  X2  and  S2  are  compared.  Yet,  in  spite  of  the 
broad  range  of  values  of  Ri,  there  is  a  tendency  (below  the 
surface  mixed  layer)  for  values  of  Logjo  (Ri)  to  remain  in 
the  neighborhood  of  zero  (i.e.  values  of  Ri  are  of  order  1). 
This  tendency  is  also  illustrated  by  the  histograms  of  Ri 
shown  in  Figure  3.7.  Hence,  some  relationship  between  the 
2m  values  of  N2  and  S2  appears  to  hold  in  a  statistical 
sense.  In  Figure  3.8  we  noted  that  N2  and  S2  values 
averaged  over  stratification  regimes  exhibited  qualitatively 
similar  behavior.  This  reinforces  the  notion  that  the  two 
parameters  are  statistically  related.  In  this  section  we 
display  the  rel at lonsrup  between  N2  and  S2  in  various 
ways,  attempt  to  explain  the  observed  relationships,  and 
assess  the  potential  for  predicting  S2  from  measurements 
of  N2. 
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4.2  PROPORTIONALITY  OF  N2  AND  S2 

Figure  4.1  is  a  plot  of  Logj  q( X2  5  versus  Log]  q 
(S2).  To  the  extent  that  density  stratification  and  vertical 
shear  are  related  by  some  power  law,  the  points  on  a  log-log 
plot  should  tend  to  lie  along  a  straight  line.-  For  example, 
if  S2  varies  as  N3  as  suggested  by  WKB  scaling  arguments 
(Grabowski,  19805,  then  the  line  will  have  a  slope  of  2/3. 
If  S2  varies  as  N2,  then  the  line  will  have  a  slope  of 
one,-  implying  that  N2  and  S2  are  directly  proportional 
and  the  proportionality  constant  is  the  Richardson  number. 
On  the  log-log  plot,  variations  m  this  proportionality 
constant  correspond  to  variations  in  the  y-mtercept  of  the 
line  with  unit  slope.-  In  Figure  4.1  diagonal  lines  with 
slopes  equal  to  one  indicate  the  various  values  of  Ri .  Note 
that  zero  values  of  \’2  or  S2  (and  negative  values  of  N2) 
are  undefined  on  this  type  of  plot.  The  horizontal  and 
vertical  lines  on  the  plot  represent  the  value  of  N2  and 
S2,  respectively,  below  which  the  signal -t o-n oi se  ratio 
(see  Section  2.1)  becomes  less  than  one. 

The  distribution  of  points  in  Figure  4.1  does  not 
clearly  indicate  any  particular  power  law  relationship. 
For  any  given  value  of  \’2  there  is  a  broad  range  of 
observed  values  of  S2.  To  compare  the  power  law  lelations 
suggested  above,  cross-correlation  coefficients  were  com¬ 
puted  between  Np  and  S2,  where  P=l,2,  and  3.  This  was 
done  for  several  YVETTE  stations  including  station  8.  The 
cross-correlation  coefficients  were  computed  for  different 
stratification  regimes  and  for  a  range  of  values  of  fz.  The 
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results  did  not  indicate  a  preferred  power  law  relation¬ 
ship,  a  conclusion  consistent  with  Figure  4.1.  It  is 
noted  that  the  points  m  Figure  4.1  do  tend  to  cluster 
around  a  line  for  which  the  value  of  P.i  is  of  order  one. 
The  distribution  shown  here  is  consistent,  with  that  shown  in 
panel  (c)  of  Figure  2.3,  except  for  the  absence  of  the 
lowest  values  of  Ri  which  are  associated  with  zero  or 
negative  values  of  N2 . 

We  now  compute,  for  every  station,  the  mean  of 
tne  2m  values  of  N2  and  the  mean  of  the  corresponding 
values  of  S2  within  each  stratification  regime.  These  mean 
values,  displayed  on  a  log-log  plot  m  Figure  4.2,  suggest 
that  N2  and  S2  are  approximately  proportional  with  a 
proportionality  constant  of  Ri=l.  The  points  in  Figure  4.2 
are  grouped  in  two  clusters.  The  cluster  with  higher  values 
of  N'2  and  S2  is  from  the  ST  regime,,  and  the  cluster 
with  lower  values  is  from  the  BST  regime.  Though  all  points 
from  the  ML  regime  have  low  values  of  N2,  they  are  not 
clustered  because  of  tne  wide  variations  m  S2. 

To  check  our  visual  perception  that  N2  is 
approximately  equal  to  S2  ,  for  each  YVETTE  station  we 
averaged  the  2m  N’2  values  and  tne  corresponding  S2  values 
over  non-overlapping  30m  intervals.  The  averaging  used 
a  triangular  weighting  scheme,  4  composite  of  the  resulting 
N2  and  S2  values  is  presented  in  Figure  4.3.  The  best 
straight-1  me  fit  to  these  points  was  determined,  using 
least  squares  regression,  to  have  a  slope  of  1.02  +  0.06 
which  is  not  significantly  different  than  1.0.  Hence, 
our  visual  perception  seems  to  be  correct. 
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S'(iz  =  2n) 


X”  versus  S'  for  a  vertical  differencing  interval  of  2m 
The  point  labels  consist  of  the  appropriate  station  number 


followed  b\  a  letter  to  indicate  the  stratification  regime 
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CORRELATION  BETWEEN  N2  AND  S2 


We  have  shown  that  the  individual  2m  values  of 
N2  and  S2  exhibit  no  obvious  relationship,  but  that  verti¬ 
cally  averaged  values  or  points  in  smoothed  profiles 
appear  to  be  better  correlated  and,  moreover,  directly 
proportional  with  a  proportionality  constant  of  one.  To 
quantify  the  correlation  between  N2  and  S2  and  the 
effects  that  data  processing  has  on  this  correlation,  we 
will  now  examine  the  linear  cross-correlation  coefficient 
between  N2  and  S2  and  the  way  in  which  it  varies  with 
vertical  differencing  interval  and  with  vertical  smoothing 
length. 


We  estimate  the  linear  correlation  between  two 
parent  populations  x  and  y  by  evaluating  the  following 
expression  for  a  sample  of  observed  pairs  '(  x  i  ,  y  j, )' , 
1=1,2,.. .M, 


r  *  (1/M),  2  (*i-jO(yi-y)/c  c  v  •• 

where  x  and  f  are  sample  averages.  The  terms  and  ~v 
denote  sample  standard  deviations.  The  sample  correlation 
coefficient  r  may  have  values  in  the  range  between  -1  and 
+1.  Values  close  to  zero  indicate  no  correlation,  while  the 
values  +  1  indicate  perfect  correlation. 

In  order  to  evaluate  the  accuracy  of  a  correlation 
coefficient  estimate,  we  must  first  evaluate  the  number  of 
statistically  independent  sample  points,  F.  Each  profile 
contains  a  different  number  of  statistically  independent 
sample  points  which  depends  on  the  profile  depth  D  and  the 
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F  =  D/t  . 

Since  N2  and  S2  may  have  different  correlation  lengths, 
we  take  *  to  be  the  geometric  mean  of  the  correlation 
lengths  of  N2  and  S2; 


Here  we  define  tx  as  the  integral  of  the  auto-covariance 
function , 

tx  =  f  lyodt  . 

o 

The  upper  limit  of  tins  integral,  tQ,.  is  defined  such  that 
Rx(t0)  =  0.1  o  which  allows  us  to  make  consistent  esti¬ 
mates  of  'x. 

We  can  assess  the  accuracy  of  correlation  coeffic¬ 
ient  estimates  by  computing  confidence  levels  at  the  68Tr 
{one  standard  deviation)  level  of  significance.  We  expand 
upon  the  procedure  outlined  by  Bendat  and  Piersol  (1971). 
Let  us  define  a  random  variable 


From  Brownlee  (1965),  w  has  an  approximately  normal  distri¬ 
bution  with  a  mean  and  standard  deiiation  of 


> 
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where  r1  is  the  correlation  coefficient  of  tne  parent 
population,  and  F  is  the  number  of  independent  observations. 


We  solve  for  r  m  terms  of  w, 

2w  , 
r  =  £ _ Z-k 

2w  ,  ,  ' 

e  +1 

The  standard  deviation  confidence  limits  r+  and 
r.  are  given  by  substituting  wi  -  for  w; 

_  ,2(w  t  t) 

-  e2(\v  t  c )  +  y 


Linear  correlation  coefficients  were  computed 
for  several  YVETTE  profiles.  The  results  for  YVETTE  sta¬ 
tion  8  are  presented  in  Figure  4.4.  Correlations  for  the 
near-surface  thermocline  (135  to  250  ra)_ ,  the  regime  below 
the  thermocline  (250  to  473ra),  and  the  combination  of 
the  two  regimes  (135  to  473n)'  are  presented.  The  \'2  and 
S2  values  were  computed  using  differencing  intervals  of  2m 
and  8m.  The  resulting  profiles  of  N2  and  S2  were 
filtered  using  triangular  weighting  functions  with  half¬ 
widths  of  5m  and  20m,  Such  a  weighting  function  is  computa¬ 
tionally  simple  and  it  has  a  well  understood  response 
function.  B  y  filtering  \ 2  and  S 2  we  increase  the 
effective  correlation  lengtus  of  f.e  profiles.  Tms  tencs 
to  reduce  the  number  of  independent  observations  and  to 
degrade  the  confidence  m  the  correlation  coefficient 
estimates . 


We  see  that  tne  correlation  m  the  combined  regime 
is  better  than  either'  of  the  separate  regimes.  This  result 


4-9 


1.0  ; 


1.0 


l 

8 


r  I 
o  r 


9a  r  i 


ST  +  BST 


TRIANGLE  FILTER 
HALF-WIDTH 

i  No  filter 
f  5m  (lov-oass) 

A  20m  (lov-pass) 

i  5m  (high-pass) 

-1  20m  (high-pass) 


Linear  correlation  coefficient  uetvec-n  X”  anti  S"  for 
the  ST  regime,  the-  BST  re-giir.e,  anti  ihe  combination  of 
the  two  regimes  at  YVETTE  station  8.  Correlation  is 
plotted  versus  vertical  differencing  interval  for 
different  triangle  filter  widths.  Confidence  limits 
are  indicated  bv  the  vertical  bars 
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applies  to  all  of  the  profiles.  This  is  due  to  the  fact 
that  in  the  near-surface  thermoclme N2  and  S2  are  much 
greater  than  in  the  regime  below.  Therefore  the  combined 
regime  contains  a  wider  range  of  values  than  either  of  the 
separate  regimes ,,  and  the  resulting  correlation  is  stronger. 
However no  general  statement  can  be  made  about  the  relative 
magnitude  of  the  correlation  coefficients  for  the  near¬ 
surface  thermocline  and  for  the  regime  below. 

Smoothing  N2  and  S2  generally  improves  correla¬ 
tions.  In  fact,  correlation  increases  as  the  smoothing 
length  increases.-  This  implies  that  long  wavelength  compon¬ 
ents  of  the  N2  and  S2  profiles  correlate  better  than 
short  wavelength  components.  We  reach  the  same  conclusion 
by  applying  a  high  pass  filter  to  N2  and  S2 .  The  high 
pass  filtered  profiles  are  computed  by  subtracting  the  low 
pass  filtered  profiles  from  the  original  profiles.  The 
resulting  frequency  response  is  the  complement  of  the 
original  frequency  response  of  the  triangular  weighting 
function.  We  see  that  the  correlation  coefficients  for  the 
high-pass  filtered  profiles  are  consistently  lower  than 
those  for  the  low-pass  filtered  profiles. 

Although  in  Figure  4.4  the  correlation  coefficient 
tends  to  decrease  with  increasing  vertical  differencing 
interval,  this  is  not  a  general  result.  There  appears  to  be 
no  consistent  relationship  between  the  magnitude  of  the 
correlation  coefficient  and  the  magnitude  of  the  vertical 
differencing  interval. 
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4.4 


PREDICTION  OF  S2  FROM  N2' 


While  Figure  4.1  shows  that  it  would  not  be 
possible  to  predict  S2  over  a  given  2m  interval  from  the 
corresponding  observed  value  of  N2,  Figures  4.2  and  4.3 
suggest  that  below  the  surface  mixed  layer  one  can  predict 
with  some  confidence  the  vertical  average  of  several  S2 
values  by  assuming  it  is  equal  to  the  vertically  averaged 
N2.  Fortunately  obtaining  a  vertically  averaged  value  of 
N2  does  not  require  a  high  resolution  density  profile  as 
provided  by  a  CTD.  It  can  be  obtained  by  simply  differenc¬ 
ing  over  a  large  (i.e.  much  greater  than  2m)  vertical 
interval  as  would  be  possible  with  Nansen  cast  data.-  In 
fact,  a  first  cut  at  the  climatology  of  small  scale  (2m) 
vertical  shear  for  the  North  Atlantic  might,  be  derived  from 
the  historical  observations  of  the  density  field  contained 
in  the  data  archives.  While  the  simple  relationship  de¬ 
scribed  above  is  probably  equally  valid  in  other  oceans, 
this  remains  to  be  demonstrated.  In  any  case,  a  relation¬ 
ship  such  as  this  can ,;  at  best,  provide  only  a  first 
approximation  to  the  mean  2m  vertical  shear  that  would 
actually  be  observed  in  any  given  profile.  This  fact  is 
demonstrated  in  Figure  4.5  which  shows  profiles  of  N2 
(solid  lines)  and  S2  (+  marks)  for  several  of  the  YVETTE 
stations.  The  points  m  the  profiles  represent  non¬ 
overlapping  50  m  block  averages  of  the  individual  N2 
and  S2  values  which  were  derived  using  a  1  z  of  2m.  This 
approach  is  similar  to  that  used  by  Simpson  (1975).  While 
the  profiles  of  verti'ally  averaged  values  are  qualitatively 
similar,  occasionally  the  agreement  at  specific  depths  is 


not  very  good.  It  may  be  that  a  temporally  averaged  value 
of  S2  '( -  z  =  2m)'  at  a  specific  depth  would  be  more  in 
agreement  with  the  local  mean  X2.-  A  single  YVETTE  profile 
is,  after  all,  only  one  realization  of  a  temporally  varying 
ocean,  and  as  such,  it  may  not  be  realistic  to  expect  better 
agreement  than  that  shown  in  Figure  4.5.  Whether  temporally 
averaged  values  would  in  fact  exhibit  better  agreement 
remains  to  be  demonstrated  with  a  time  series  of  velocity 
profiles. 

Ir  Figures  3.5  and  3.6  we  showed  the  effect  on 
S2  of  increasing  the  vertical  differencing  interval. 
This  effect  is  shown  again  m  the  sequence  of  T2  versus 
S2  plots  presented  m  Figure  4.6.-  4s  expected  from  the 
previous  discussions,  decreases  and  the  ratio  of  X2 
to  S2  increases  as  Lz  increases  (i.e.  the  relationship 
between  N2  and  S2  changes  as  Lz  changes).  However, 
i  Figure  4.6  also  reveals  that  the  scatter  of  the  points  in 

each  cluster  increases  as  Lz  increases.  Hence,-  for  a  given 
value  of  X2  we  observe  an  increasing  percentage  variation  of 
S2.  This  suggests  that  we  can  predict  a  mean  value  of 
>  S2  over  a  2m  vertical  scale  with  a  smaller  percent  uncer¬ 

tainty  that  we  can  predict  a  mean  value  over  a  larger  verti¬ 
cal  scale.  The  reason  for  this  is  not  yet  clear.  Nonethe¬ 
less,  the  empirical  relationship  between  N?  and  is  sensi- 
1  live  to  the  magnitude  of  Lz  and  we  have  more  confidence  in 

the  relationship  if  Lz  is  small.  However,  we  noted  in 
Figure  4.5  that  even  for  _z  =  2m,  the  relationship  is  only 
approximate. 

) 
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In  Figure  3.8  we  showed  that  S2  for  each  stra¬ 
tification _ regime  varied  in  almost  parallel_  fashion  with 

N2.  If  N'2  were  exactly  proportional  to  S2,,  then  their 
ratio  would  be  constant.  In  Figure  4.7  we  present  a  plot  of 
the  ratio  of  N2  to  S2  versus  station  number.  The  stations 
are  arranged  according  to  decreasing  S2  as  in  Figure  3.8. 
The  ratio  for  each  stratification  regime  is  plotted  as  well 
as  the  ratio  for  the  whole  profile.  To  the  extent  that  the 
ratio  is  not  constant,  the  proportionality  relationship 
fails.  We  note  that  while  the  ratio  below  the  surface  mixed 
layer  is  usually  of  order  one,  the  actual  value  varies  from 
station  to  station.  The  high  ratio  observed  at  stations  F6 
and  N6  results  from  anomalously  low  values  of  S2.  The 
significance  of  this  is  not  yet  clear,-  but  it  suggests  that 
the  relationship  between  \’2  and  $2  which  holds  in  the 
western  North  Atlantic  may  "ot  be  universal.  There  is  a 
tendency  seen  in  Figure  4.7  for  different  stratification 
regimes  to  vary  in  parallel  fashion,  that  is,  if  the  ratio 
is  relatively  high  in  the  ST  regime,  it  is  usually  high  in 
the  BST  regime  also.  Hence,  whatever  it  is  that  causes  the 
ratio  to  vary  from  one  station  to  another,  it  tends  to 
affect  the  whole  profile  in  the  same  way.  As  noted  earlier 
a  temporally  averaged  S2  at  each  station  might  exhibit 
better  agreement  with  the  local  \'2. 

The  ratio  of  N2  to  S2  in  Figure  4.7  is  almost 
always  higher  in  the  ST  . eg tme  than  below.  This  might 
possibly  be  explained  in  terms  of  the  natural  variation  in 
vertical  scales  that  occurs  in  the  stratified  ocean.  An 
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Figure  4.7  Ratio  of  N*  to  S"  versus  station  number  for  the 
different  stratification  regimes. 
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internal  wave  propagating  vertically  through  the  water 
column  decreases  its  vertical  wavelength  as  N2  increases. 
41so  the  Osmtdov  scale  of  turbulent  overturning  decreases  as 
\'2  increases.  Hence,  within  the  ST  regime,  variations  in 
horizontal  velocity  tend  to  occui  over  smaller  vertical 
scales  than  in  the  BST  regime.  With  any  fixed  vertical 
differencing  interval,  such  as  £z  =  2m,  we  are  not  able  to 
resolve  variations  that  occur  over  smaller  vertical  scales. 
Hence,  our  estimate  of  s"^  will  always  be  somewhat  too 
low,  but  it  will  be  proportionately  lower  within  the 
ST  regime  than  below.  This  is  consistent  with  the  observa¬ 
tion  that  the  ratio  of  "ic2  to  is  usually  higher  within 
the  ST  regime. 

In  Figure  4.7  there  appears  to  be  a  slight  trend 
of  increasing  ratio  from  left  to  right  across  the  figure. 
Recall  that  both  N2  and  S2  are  generally  decreasing  in 
this  direction,  and  hence  the  increasing  ratio  suggests  that 
S2  must  be  decreasing  faster  than  N'2.  To  the  extent  that 
the  perceived  trend  is  significant,  the  simple  proportional¬ 
ity  relationship  between  \2  and  S2  fails  again. 

Intuitively  it  is  perhaps  more  surprising  that 
N 2  and  S2  are  as  consistently  related  to  each  other 
from  station  to  station  as  the  observations  indicate,  than 
the  fact  that  the  relationship  is  not  perfect.  This  is 
hpcause  it  seems  reasonable  to  expect  that  the  S2  level 
within  any  given  profile  is  dependent  to  some  extent  on  the 
variable  influence  of  present  and  past  forcing  events.  In 
other  words,  two  stations  with  the  same  W2  levels  might 
have  different  S2  levels  because  of  a  difference  in  recent 
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or  local  forcing  conditions.  As  the  kinetic  energy  is  being 
put  into  the  water  column,  it  will  tend  to  arrange  itself  in 
accordance  with  the  local  density  profile,  e.g.  vertical 
shear  will  be  greatest  in  regions  of  high  density  gradient. 
In  the  case  of  YVETTE  stations  8  and  9,  we  note  that  they 
were  occupied  at  the  same  location  and  were  separated  in 
time  by  approximately  10  hours  (Table  1.1)’,-  The  N  2 
profiles  for  the  two  stations  are  almost  identical,  but  the 
S2  profiles  are  different  (cf.  Figures  2.3  and  A. 8),.  with 
YVETTE  8  having  a  higher  S2.  This  accounts  for  the 
lower  ratio  of  N 2  to  S2  exhibited  by  YVETTE  8  than 
by  YVETTE  9  in  Figure  4.7.  Hence  the  N'2  value  for  a 
particular  portion  of  the  water  column  may  establish,  in  a 
statistical  sense,,  an  "expected"  value  for  S2.-  but  whether 
or  not  the  observed  S2  level  fulfills  these  expectations 
depends  on  the  availability  of  sources  of  kinetic  energy. 
The  trend  of  increasing  ratio  seen  in  Figure  4.7  might  be  an 
artifact  of  arranging  the  stations  in  order  of  decreasing  S2 
which  automatically  incorporates  the  effects  of  decreasing 
availability  of  kinetic  energy. 
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Section  5 

SUMMARY  AND  CONCLUSIONS 


Vertical  profiles  of  density  and  horizontal 
velocity  from  a  diverse  range  of  oceanic  conditions  have 
been  examined.  Vertical  gradients  of  these  quantities  are 
computed  by  taking  simple  differences  over  a  vertical 
interval,  iz.  The  magnitudes  of  the  vertical  gradients, 
in  particular  Brunt-Vai  s'al’a  frequency  squared  (N2)  and 
vertical  shear  squared  (S2),  appear  to  be  arranged 
into  distinct  regimes  of  differing  stratification.  The 
surface  mixed  layer  (ML),  sampled  in  some  profiles,  is 
characterized  by  vertically  homogeneous  density  and  weak 
mean  vertical  shear.  The  near-surface  thermorline  (ST)  is 
characterized  by  a  broad  range  of  values  in  both  \’2  and 
S2.  The  regime  below  the  ST  (the  BST  regime)  is 
characterized  by  uniformly  small  values  of  N2  and  S2. 
Vertical  profiles  of  Logjo  (Ri),;  where  Ri  *  N2/S2,, 
exhibit  values  spanninc  several  orders  of  magnitude, 
but  show  only  very  slight  differences  between  the  ST  and  BST 
regimes.  The  ML,  on  the  other  hand,  exhibits  very  low  Ri 
values  and  is  therefore  distinctly  different  than  the 
deeper  regimes.  The  lowest  values  of  Ri  generally  occur 
where  N2  is  low  rather  than  where  S2  is  high. 

Histograms  for  eat  h  stratification  regime  reteal 
that  the  mean,  mode,  and  variance  of  S2  distributions 
decrease  with  increasing  ~z.  This  is  because  the  highest 
values  of  S2  are  associated  with  small  scale  fluctua¬ 
tions  in  the  velocity  profile.  As  Lz  increases,  the  ability 
to  resolve  these  small  scale  fluctuations  is  lost.  The 
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decrease  in  mean  S  2  (  S  2  }  with  increasing  :z  is  con¬ 
sistent  enough  to  suggest  the  possiblity  of  extrapolating  to 
vertical  scales  which  are  not  well  resolved  by  YVETTE.  This 
might  allow  direct  intercomparison  of  S2  levels  measured 
by  sensors  with  different  vertical  resolution  capabilities.- 

Since  the  S2  distributions  are  sensitive  to 
changes  in  £z,  the  same  is  true  for  distributions  of  Ri . 
The  Ri  values  become  generally  larger  as  £ z  increases.- 
While  the  S2  distributions  for  the  ST  and  BST  regimes  are 
distinctly  different,  the  corresponding  Ri  d 1 st r iubi tons  are 
very  similar.  Ri  distributions  are  more  sensitive  to 
changes  in  iz  than  to  changes  in  stratification  regime. 
This  suggests  that  the  details  of  the  shape  of  the  Ri 
distributions  may  depend  more  on  the  vertical  resoluton 
capabilities  of  the  velocity  sensor  than  on  the  ocean 
environment . 

In  the  North  Atlantic  iT2  and  S2  appear  to 
be  approximately  proportional.  if  iz  =  2m,  the  propor¬ 
tionality  constant  is  Ri  =  1.  Both  N 2  and  S2  tend  to 
decrease  with  distance  from  the  Gulf  Stream.-  Both  appear  to 
increase  toward  the  center  of  cyclonic  mesoscale  eddies.-  A 
first  approximation  to  the  climatology  of  small  scale 
vertical  shear  in  the  North  Atlantic  might  be  derived 
from  the  historical  observations  of  the  density  field 
contained  m  the  data  archives.  The  relationship  between 
N2  and  S2  in  other  ocean  areas  needs  to  be  established 
with  additional  YVETTE-type  observations. 

The  correlation  coefficient  between  N2  and  S2 
values  (computed  over  a  specified  ,‘z)‘  increases  with  in¬ 
creasing  vertical  smoothing  interval  (vertical  filter 
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length)'.  Hence  the  correlation  between  long  vertical 
wavelength  components  of  the  profiles  is  better  than  the 
correlation  between  the  short  vertical  wavelength  compo¬ 
nents.  Apparently,  small  scale  vertical  shear  has  an 
expected  value  which  is  modulated  by  the  large  scale  density 
gradient.  Since  the  magnitude  of  S2  depends  on  vertical 
changes  in  either  the  magnitude  or  direction  of  the  hori¬ 
zontal  velocity  vector,  the  modulation  by  the  large  scale 
density  gradient  must  be  associated  with  a  modulation  of 
either  nr  both  of  these  quantities.  Changes  in  speed 
correspond  to  changes  in  the  kinetic  energy  level:  changes 
in  direction  might  be  due  to  internal  waves.  The  YVETTE 
profiles  show  that  changes  in  both  speed  and  direction 
contribute  to  changes  m  S2.  The  correlation  coefficient 
exhibits  no  consistent  trend  with  changes  in  the  vertical 
differencing  interval. 

The  proportionality  constant  between  N2  and 
S2  increases  as  iiz  increases.  This  means  that  the 
relationship  used  to  predict  S2  from  measured  values 
of  N2  depends  on  the  particular  vertical  differencing 
interval  employed..  Our  analysis  of  the  YVETTE  data  suggests 
that  S2  over  a  2m  interval  can  be  predicted  with  a  smaller 
percent  uncertainty  than  S2  over  a  larger  vertical 
i  n  terval  .• 


Although  a  proportionality  relationship  between 
N2  and  S2  seems  to  be  a  good  first  approximation, 
S2  is  clearly  dependent  to  some  extent  on  other  factors, 
such  as  the  availability  of  sources  of  kinetic  energy.  It 
is  possible  that  temporally  averaged  values  of  S2  would 


exhibit  better  agreement  with  the  local  N2  profile.  The 
apparent  difference  between  the  ratio  of  N2  to  S2  for 
the  ST  regime  and  the  ratio  for  the  BST  regime  may  be 
due  to  the  natural  variation  of  vertical  scales  (of  internal 
waves  and  the  Osmidov  length)  that  occurs  in  the  stratified 


ocean. 


Section  6 

RECOMMENDATION’S  FOR  FURTHER  WORK 


Although  a  diverse  range  of  dynamically  differenr 
oceanic  conditions  has  been  examined  in  this  data  set,  most 
of  the  stations  are  located  in  the  western  North  Atlantic. 
The  anomalously  low  levels  of  small  scale  shear  observed  in 
the  Norwegian  Coastal  Current  and  the  equatorial  Pacific 
cause  these  stations  to  deviate  most  from  the  more  or  less 
consistent  pattern  established  by  the  other  stations.  To 
determine  the  generality  of  the  conclusions  of  this  report.- 
other  ocean  areas  (e.g.  North  Pacific  and  Indian  Ocean)  and 
mesoscale  features  (e.g.  anticyclonic  eddies  and  fronts) 
need  to  be  sampled  with  an  YVETTE-type  profiler. 

Each  profile  examined  in  this  report  provides  a 
single  realization  of  a  temporally  varying  ocean.-  It  is  to 
be  expected  that  the  levels  of  small  scale  shear  at  any 
location  will  vary  in  time  as  a  result  of  forcing  events, 
inertial  oscillations,  etc.  Therefore,  we  also  need  to 
obtain  and  examine  time  series  of  profiles  in  several 
locations  to  (1)  determine  how  temporal  variability  affects 
the  relationship  between  N2  and  $2  and  (2)  provide 
enough  data  to  resolve  the  vertical  wavenumber  spectrum  for 
the  different  stratification  regimes  and  different  ocean 
areas . 


A  climatology  of  large  scale  vertical  density 
gradients  could  be  derived  from  the  historical  observa¬ 
tions  of  the  density  field  which  are  contained  in  the  data 
archives.  From  this  characterization  a  first  cut  at  a 
climatology  of  small  scale  shear  might  be  derived. 
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Appendix  A 
YVETTE  PROFILES 
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Eigure  A.l  Profiles  of  (a)'  temperature,  (b)‘  salinity,  (c)  density 

(sigma-t),  (d)  Brunt-Vaisala  frequency,  and  ,(e)  vertical 
gradient  of  temperature  lor  YVETTE  station  5  occupied  on 
5  November  1975  at  32°19'N.  64°34'\v  (near  Bermuda).  The 
profiles  can  be  divided  into  stratification  regimes  as 
indicated  by  the  dotted  lines  at  55  and  150  dbars 


Figure  A. 6  As  m  Figure  A. 3  except  for  YtETTE  station  9 
m  the  Sargasso  Sea.  The  dotted  lines  are  at 
130  and  280  dbars. 


Figure  A. 12 


As  in  Figure  A. 3  except  for  YVETTE  station  11 
in  the  Gulf  Stream.  The  dotted  line  is  at 
150  dbars. 


Figure  A. 13  As  m  Figure  A.l  except  for  YVETTE  station  12 

occupied  on  10  November  1975  at  35'15'N,  69C07'K 
(Gulf  Stream).  The  dotted  line  is  at  135  dbars. 


Figure  A. 14  As  in  Figure  A. 2  except  for  YVETTE  station  12 
in  the  Gulf  Stream.  The  dotted  line  is  at 
135  dbars. 


Figure  A. 17  As  in  Figure  A. 2  except  for  YVETTE  station  18 
near  the  Bahamas.  The  clotted  line  is  at 
225  dbars. 


Figure  A. 18  As  in  Figure  A. 3  except  for  YVETTE  station  18 

near  the  Bahamas.  The  dotted  line  is  at  225  dbars 


Figure  A. 19  As  m  Figure  A.l  except  for  YVETTE  station  21 
occupied  on  9  May  1977  at  22'I27’N,  70’37'Vf  (a 
center  of  thermocline  eddy  near  the  Bahamas). 
The  dotted  line  is  at  280  dbars. 


Figure  A. 23  As  in  Figure  A. 2  except  for  YVETTE  station  23  on 
the  outer  edge  of  a  Gulf  Stream  ring.  The  dotted 
line  is  at  120  dbars. 


» 


Figure  A. 24  As  m  Figure  A. 3  excpet  for  YVETTE  station  23 
at  the  outer  edge  of  a  Gulf  Strean  ring.  The 
dotted  line  is  at  120  dbars. 
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Figure  A. 26  As  in  Figure  A  2  except  for  YVETTE  station  24 
midway  along  tne  radius  of  a  Gulf  Stream  ring. 
The  dotted  line  is  at  82  dbars. 
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Figure  A. 34 


As  m  Figure  A.l  except  for  YVETTE 
occupied  during  1973  in  a  .Norwegian 
dotted  line  is  at  30  dbars. 
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Figure  A. 40  As  in  Figure  A.l  except  for  YVETTE  station  E6 
occupied  in  July  1979  at  0C03'N,  109°57Mf 
(equatorial  Pacific).  The  dotted  lines  are  at 
115  and  240  dbars.  In  this,  case  the  upper  regime 
is  referred  to  as  the  near-surface  thermoclme  (Sv) 
the  central  regime  is  a  subsurface  pycnostad  (P). 
and  the  lower  regime  is  below  the  pycnostad  (BP). 


'igure  A.  2  except 
equatorial  Pacifi 
115  and  240  dbars 


Appendix  B 
YVETTE  HISTOGRAMS 
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ST  (21m  -  160m)  .3  Mean  -  1.803 

No.  Obs.  -  142  Median  =  1.153 

Moan  -  1.862  Variance  =  3.508 

Median  -  1 . 225 
Variance  --  1.035 
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